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Abstract. 

We have developed an excitation model for small molecules 
including radiative pumping by dust photons from near in¬ 
frared to submillimeter wavelengths. This model applies to 
molecules within bright photodissociation regions in galactic 
star-forming regions and starburst galaxies. In such environ¬ 
ments, the near infrared photons emitted by polycyclic aro¬ 
matic hydrocarbon molecules (PAHs) or small carbon grains, 
5 /rm < A < 20 /rm, are able to penetrate molecular regions 
and pump the molecules in an excited vibration state. The 
far infrared and submillimeter transitions involved in the de¬ 
excitation cascade are strongly affected by this process. Their 
intensities can be enhanced by several orders of magnitude. We 
have applied this model to H 2 O and NH 3 , ortho and para 
species. The behaviour of the FIR rotation lines with respect 
to the gas density and the molecule column density is strongly 
modified compared to pure collisional models. At moderate 
densities, nH 2 < 10® cm~^, radiative pumping dominates the 
excitation of the rotation ladder, and the FIR lines are good 
probes of the molecule column density. The near infrared ab¬ 
sorption and emission lines predicted to appear between 6 and 
7 fim for H 2 O and from 10 to 12 /im for NH 3 can also be 
used for this purpose. Concerning ortho — H 2 O, a quasi reso¬ 
nant pumping of the V 2 6.18/rm transition by photons of the 
6.2 nm PAH band occurs. A strong de-excitation emission line 
is expected at 6.64 fim. 

Key words: dust - ISM: molecules - ISM:clouds - reflection 
nebulae - Infrared:ISM:lines and bands - Galaxiesistarburst 


1. Introduction 

Modelling the formation of molecular rotation lines in inter¬ 
stellar clouds is a crucial task both for theoretical and obser¬ 
vational astrophysics. On t he theoretical side, early stud ies by 
Goldreich and Kwan (1974) and Goldsmith and Langer (1978) 
have shown that the rotation lines from polar and/or abun¬ 
dant molecules like CO, O 2 , HD, H 2 O and other hydrides 
should be the major cooling agents of the molecular gas. To¬ 
gether with the heating mechanisms (e.g. cosmis rays, H 2 for- 
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mation, gravitational contraction and grain heating), they con¬ 
trol the thermal balance of dense clouds, and thus, the condi¬ 
tions that will allow the contraction toward star formation. On 
the observational side, millimeter and submillimeter molecular 
lines are used to probe several physical parameters of molec¬ 
ular clouds, such as the gas temperature, the H 2 density, the 
abundance of molecules, and the gas dynamics. Densities and 
temperatures in molecular clouds are such that the collision 
rate with H 2 molecules generally does not allow to reach the 
thermodynami c equ ilibrium of the rotation ladder (see Goldre¬ 
ich and Kwan 1974 ). Instead, the coupling with the radiation 
field is important, both via a bsor ption and spontaneous emis¬ 
sion. Takahashi et al. (1983 and 1985), have shown how far 
infrared dust photons, A = 30 to 100/rm are able to pump 
the rotation ladder of the H 2 O m olecul e within warm molecu¬ 
lar cores. Carroll and Goldsmith (1981) have investigated the 
pumping of molecules via their near infrared (NIR) vibration 
transitions. At that time it was assumed that a significant frac¬ 
tion of the dust emission could be in the NIR range, only in the 
immediate vicinity of hot stars. However, the effect of the NIR 
radiation attributed to large polycyclic aromatic hydrocarbon 
molecules (PAHs) and very small grains, wh ich m akes up to 
30 % of the total dust emission (Puget et al. 1985), has never 
been taken into account. This radiation is characterised by a 
strong continuum and broad emission features, located at 6.2 
and 7.7 pm for the strongest ones, wavelengths which coincide 
with the vibration transitions of some molecules. For instance 
the v2 transition of H 2 O and the w4 transition of NH 3 occur 
at the maximum of the 6.2 pm PAH feature, and the 7.8 pm 
vibration of GS falls within the wavelength range of the 7.7 
pm feature. There is now clear observational evidence that the 
NIR emission bands are not only present in peculiar reflexion 
nebulae, but are a general component of the dust radiation: 
detection of the 3.3 and 6.2 pm ba nds in our galaxy by Gi¬ 
ard et al. 


1994a and Ristorcelli et al. 


1994 


and, more recently, 
the m easur ements perform ed by the ISO satellite (Bo ulange r 
et al. 1996| , Mattila et al. 1996| , and Verstraete et al. 199(| ). 
The aim of this paper is to revisit the problem of the coupling 
between molecules and the dust radiation field, in the frame¬ 
work of a model which includes the PAH emission bands. The 
mechanism, and the rates involved in some representative as- 
trophysical environments, are presented in Sect. A detailed 
model for the molecular excitation in star-forming photodisso¬ 
ciation regions (PDR) is develloped in Sect. ^ The results for 
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two molecules, H 2 O and NH 3 , are presented in Sect. 


2 . Near infrared pumping of rotation transitions 

Near infrared (NIR) PAH photons, continuum or bands, can 
excite the vibrations of a molecule. This absorption does also 
imply some rotational excitation (AJ 7 ^ 0). After spontaneous 
de-excitation of the vibration mode, which occurs almost im¬ 
mediately {Aij ~ 10s~^), the molecule can be left in an ex¬ 
cited rotation state. Subsequent excitation from this state to 
an upper one can follow via the same mechanism. We have 
listed in Table some molecules of astrophysical interest which 
may be affected by this process, because one of the vibration 
transitions happen to be at a wavelength larger than 6 /rm 
(< 1670cm“^), where the PAH emission is strong. 


Table 1 . Some molecules of astrophysical interest which hap¬ 
pen to have one or several vibration transitions that can be 
pumped by PAH photons. Frequencies and Einstei n coe ffi- 
cients are from the compilation of Encrenaz et al. 1992. If 
not available the Einstein coefficient has been approximated 
to 10s“^. Column 4: indicative abundance relative to hy¬ 
drogen. When available, we have used the value fr om t he 
dust/gas chemical model of Shalabiea and Greenberg 1994 at 
a time evolutio n of 1 0^ years. The abundance of is from 

Phillips et al. 


1993 


Column 5: optical depth in the absorb¬ 
ing vibration transition for a typical PDR column density: 

2Nh2 = 10 ^*^ 


Molecule 

^vib 

cm~^ 

-^vib 

S~^ 

Abundance 

T 

H2O 

U 2 1620 

20 


19 

NH3 

1^2 950 

16 

10 "® 

7 


Ui 1627 

5 


0.5 

CO2 

1^2 670 

1.4 

10 "® 

2 

C3H2 (cyclic) 

j /3 1277 

10 

2 10 “® 

0.04 

CH4 

Vi 1306 

2.6 

10 "® 

0.5 

HCN 

1^2 712 

1.7 

10 "® 

2 

HC3N 

123 663 

4.7 

10 "® 

0.006 

H3O+ 

1^4 1626 

10 

3 10“® 

0.003 

CH3OH 

1 / 4 ...Us, Vio, V\\ 

10 

3 10“^ 

1.2 


1033 ... 1478 




CS 

1285 

10 

3 10“® 

5.5 


where (respectively gi, Ei) are the statistical 

weight and the energy of the upper (respectively lower) level. 

The radiative excitation rate is given by the product of the 
Einstein coefficient by the spectral radiation density: 


^ ~ 6.8 10 "^^ 


< E> , 

MJy/sr 


V 


-3 


( 3 ) 


where A^i is the spontaneous emission rate, < E > is 
the average spectral density of the radiation field in the line 
(average over frequency and solid angle). The spectral den¬ 
sity at the maximum of the 6.2 and 7.7 gm PAH features 
is a few 10“^ MJy/sr within typical photodissociation regions 
(PDRs) a ssocia ted t o star forming complexes (e.g. M17, see Gi- 


ard et al. 


1992 


and 


19941 ,) or starburst galaxies. This implies 


that in such regions and their vicinity, the radiative excita¬ 
tion rate to the vibration excited state, ~ 3 10“® s“^, can be 
higher than the collision excitation rate of the rotation ladder, 
~ 7 10"^ Table I compares the radiative and collisional 
rates obtained in the case of the excitation of the 221 rotation 
level of ortho — H 2 O, 111 cm~^ above the ground state. Fig. 
1 ^ shows in detail how this level can be excited via radiative 
excitation at 6.18 fim of the v2 lio state, and spontaneous 
emission, Aij ~ 10We have considered different astro- 
physical environnements where radiation can play a significant 
role: a typical ^ 1/772 galactic PDR in a regio n of massive 
star formation (see Tielens and Hollenbach 1985 ); a molecular 
cloud nearby to this PDR (e.g. M17 SW); the central regions 
of a starbust galaxy (e.g. M82); and cold molecular gas within 
the galactic molecular ring. In the last case, the low value of 
the radiation field is balanced by the poor collisional efficiency 
at low temp eratur es. T = 50 K. The collision rates are from 
Green et al. 1993 . The figures in this table are only given to 
illustrate the fact that in such astrophysical environments ra¬ 
diative excitation by NIR PAH photons efficiently competes 
with collisional excitation. A correct treatment has to take into 
account the escape probability of the photons from the cloud. 
This will be done in the next section. 


Table 2 . Excitation of the 221 rotation level of ortho — H2O, 
111 cm~^ above the ground level, for different astrophysical 
environment s, fcj) , direct collisional rate from the ground level 
(Green et al. 1993). k~^^, radiative excitation rate via the v 2 lio 
vibration state (absorption of a 6.18/rm photon). 


In order to get a first quantitative estimate of the compe¬ 
tition between the proposed mechanism (radiative NIR pump¬ 
ing) and collisional excitation, we have to figure out what are 
the rates of each elementary process. Concerning the effect of 
collisions with H 2 , t he dee xcitation rate is approximately (see 


Goldreich and Kwan 


1974) : 



< E > 
MJy/sr 

T 

K 

riH^ 

cm~^ 

Ko\ 

S ^ 

^rad 

S~^ 

PDR 

TIF 

120 


6.4 10 ”® 

1.8 10 "® 

M17 SW 

1.5 10“* 

50 

10® 

3.7 10 “® 

0.9 10 "® 

Starburst 

10'‘ 

120 

10“^ 

6.4 10“® 

0 

1—1 

0 

1 

■SI 

Gal. 772 

40 

50 

10“* 

3.7 10“^° 

2.4 10"® 


UHn 


r.—1 7^ 

^col.deexc. ^ ^ ^ \ 0-5 _ 

s~^ 1007s' 10^ cm“® 


( 1 ) 


where T is the gas temperature, and njfj the gas density. 
The excitation rate follows from the detailed balance: 


^ k-o\.a...c.^exp{-{E^ - Ei)/kT) 
s yi 


( 2 ) 


3 . The model 

Our model solves the equations of statistical equilibrium for a 
given molecule, including rotation levels in the ground and first 
excited vibration states, assuming an homogenous gas cloud (ie 
the excitation of the molecule doesn’t vary across the cloud) : 
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71^2 ^ij') ^ ^ ^j'i-') 

0 (4) 

where Xi is the molecular fraction in each energy level 
(^^Xi = 1), Cij (cm~^s~^) are the collision rates, Rij (s~^) 
the radiative rates and the photon frequency of the transi¬ 
tion. The radiative transition rates read: 


dxi 

dt 


Rij = Aij{l + < R >) for i > j 

c2 

Rij = Aji ——— — < h > for i < j 


‘ 2huf^ Qi 


( 5 ) 




<T>: Average line optical depth 

Fig. 1. Escape probabilities for molecule and dust photons versus 
the line average optical depth. 


< Iv > IS the averaged spectral density of the radiation 
held at the frequency of the transition. This average has to be 
performed both over the velocity distribution of the molecules, 
and over the solid angles. 

Concerning the radiative transfer, we assume that the PDR 
is a clumpy medium at small scale, with molecular clumps 
embedded in a lower density interclump gas, hh = 10® to 
10^ cm~^. There are several observational evidences which sup¬ 
port this hypothesis, among which one can cite: 1/ the high Cl 
and CII column densities measured to ward the M17-SW PDR 
and molecular cloud by Genzel et al. 


1988 


and Stutzki et al. 
1986, and 2/ the extension of the PAH emission in the Orion 


bar, which shows that the UV flux penetrates much deeper into 
the molecular cloud than what is allow ed by extinction in a ho - 
mogeneous gas (see Tielens et al.l993 and Giard et al. 1994b). 
We assume in the following that the UV radiation is converted 
to near infrared PAH photons in the interclump medium and 
the outer parts of the molecular clumps, the molecules be¬ 
ing at larger Av behind the dissociation limit. The peculiarity 



infrared exciting radiation (the PAHs) is not at thermal equi¬ 
librium. Here we didn’t try to take into account the absorption 
and emission by dust in a selfconsistent way. Instead, we have 
neglected the absorption by dust an d we u se for the dust radi¬ 
ation field the model of Desert et al. 1990, with the flu x scaled 
on the 3.3 jim measurements of Giard et al. 


1994b on M17 


South-West and the Orion bar. The only extinction is then the 
one of the molecules into the clumps. We can use for the dust 
photons a radiative transfer of the form {l — exp{—r))/T, which 
has been shown to be a good approximation for a medium 
with absorbing clumps and t ransp arent intercloud phase (see 


spai 

Appendix A in Giard et al. 1992| ). However, the uncertainty 


introduced in the model by the transfer used for the dust pho¬ 
tons is attenuated by the integration over the line prohle which 
follows in Equ. ^ below. 

With such hypothesis the trausfer equation reads: 


1 

91^-u. 


( 7 ) 


u and 1 denote the upper and lower level of the transition, 
and T is the wavelength dependent molecular optical depth: 

t{v) = Nxu<l>{v)^^ ( 8 ) 

N is the molecular column density, and the normalized 
doppler profile resulting from the velocity distribution of the 
molecules. In the most general case of our aproximations, the 
profile averaged radiation density used in Equ. (statistical 
equilibrium) reads: 


< Iv >— 5(1 — Pm) H— 


( 9 ) 


where the escape probabilities for molecule and dust pho¬ 
tons read: 


Pm{<T>) = J exp{-T{v)/2)(t>{v)dv 

M< r». J (10) 

The factor 2 in Equ. ^ and 0 arises from the fact that, on 
average, the column density surrounding any point within the 
cloud is half that seen by the observer. 

Pm and Pd have been computed numerically for a Maxellian 
velocity distribution. They are plotted in Fig. as a function 
of the average line optical depth: 


N c 


Aui 

47r 


hu 

'J 


( 11 ) 


where DV is the width of the projected velocity distribu¬ 
tion and S the source function. 

Finally, the integrated intensity of the lines are computed 
by evaluating the power radiated by the molecules. Normalized 
to the velocity width it reads: 


R = 5(1 - exp(-T)) -b Fd (6) 

r 

where Fd is the dust and PAHs spectral density in the 
continuum, and S is the line source function: 


_ FR _ 

Wm~^s'i — '^(km/s)~'^ 


- < r > {PmS - ^Pd) 
c 2i 


(12) 
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Table 3 . Calculated population fractions of the lower ro¬ 
tational levels of H 2 O (T = 120 K , uh^ = 10 “* cm“®, 

NH2O = 10^® cm~^{km/ s)~^, ortho/para = 3). Colnmn 4: ra¬ 
tio of the level population fraction to the population obtained 
with a pure collisional model. 


a) ortho — H2O 


Energy 

state 



fraction 

fraction 

[ cm~^] 

V 

JK 

-K+ 

in state 

/fraction(Coll.) 

23.8 

000 

1 

0 

1 

8.75e-01 

1.0 

42.4 

000 

1 

1 

0 

9.18e-02 

0.9 

79.5 

000 

2 

1 

2 

3.21e-02 

2.5 

134.9 

000 

2 

2 

1 

2.59e-04 

100.3 

136.8 

000 

3 

0 

3 

3.58e-04 

25.9 

173.4 

000 

3 

1 

2 

1.53e-05 

66.4 

212.2 

000 

3 

2 

1 

1.20e-05 

225.5 

224.8 

000 

4 

1 

4 

1.19e-07 

4.8 

285.4 

000 

3 

3 

0 

4.47e-08 

656.7 

b) para 

-H 20 





Energy 

state 



fraction 

fraction 

[ cm~^\ 

V 

JK 

-K+ 

in state 

/fraction(coll.) 

0.0 

000 

0 

0 

0 

8.68E-01 

0.9 

37.1 

000 

1 

1 

1 

1.19e-01 

2.0 

70.1 

000 

2 

0 

2 

1.18e-02 

15.7 

95.2 

000 

2 

1 

1 

6.35e-04 

237.6 

136.2 

000 

2 

2 

0 

5.63e-05 

136.9 

142.3 

000 

3 

1 

3 

1.08e-05 

37.7 

206.3 

000 

3 

2 

2 

1.35e-07 

251.4 

222.1 

000 

4 

0 

4 

2.55e-08 

1.2 

275.5 

000 

4 

1 

3 

1.12e-09 

2.7 

285.2 

000 

3 

3 

1 

2.63e-08 

648.0 


We have run the model for H2O and NH3, two polar 
molecules likely to be abundant in PDRs because they are ex¬ 
pected to be formed in grain mantles and released to th e gas 
in such regions (see e.g. Shalabiea and Greenberg 1994). The 
physical parameters of the regions are typical of a PDR in a 
region of high mass star formation (first line in Table |^. Con- 
cernin g the dust spectrum, we have used the model of Desert 
et al. (1990) in the case of dust irradiation by an O star. The 
dust flux has been scaled so that the peak 6.2 /rm brightness is 
3 10^ MJy/sr. This is the value expected for this type of PDRs, 
on the basis of the 3.3 feature flux meas ured to ward M17 
south-west and the Orion bar (G iard et al. 1994b|) .The ISO- 
SWS spectrum of Verstraete et al (1991) on the interface region 
in M17 is about a factor of 15 below the above surface bright¬ 
ness, because of dilution in the ISO-SWS aperture (14" x 20") 
and averaging over 3 different positions. 

The coupled system of Equs. and statistical equilib¬ 
rium and radiative transfer, have been solved numerically for 
each molecule, varying the hydrogen density and the molecule 
column density. We selected the 16 lower rotation levels in 
each vibration state of the molecule studied. The collision de¬ 
excitation rates for transitions between rotation lev els in the 
ground vibration state a re tak en from Green et al. (1993) for 
H2O, and Danby et al. (1988) for NH3. Concerning the col¬ 
lision deexcitation rates between different vibration states, we 


Table 4 . Calculated population fractions of the lower ro¬ 
tational levels of NH 3 (T = 120K, nH 2 = 10“* cm“®, 

N NH2 = 10^® cm~^{km/ s)~^, ortho/para = 2). Same as Table 

I 

a) ortho — NH3 


Energy 

state 

fraction 

fraction 

[ cm“^] 

V JK± 

in state 

/fraction (coll.) 

0.4 

0000 0 0- 

2.59e-01 

0.9 

19.5 

0000 1 0+ 

7.95e-02 

1.8 

60.0 

0000 2 0- 

1.19e-03 

10.6 

85.9 

0000 3 3- 

3.04e-01 

1.0 

86.7 

0000 3 3-b 

3.07e-01 

1.0 

118.8 

0000 3 0-b 

1.57e-06 

3.8 

165.3 

0000 4 3- 

4.09e-04 

11.2 

166.1 

0000 4 3-b 

3.89e-04 

18.0 

b) para 

-NH3 



Energy 

state 

fraction 

fraction 

[ cm~^] 

V JK± 

in state 

/fraction(coll.) 

16.2 

0000 1 1- 

2.65e-01 

1.0 

17.0 

0000 1 1+ 

2.45e-01 

1.0 

44.8 

0000 2 2+ 

1.97e-01 

1.0 

45.6 

0000 2 2- 

1.80e-01 

1.0 

55.9 

0000 2 1- 

9.80e-04 

3.2 

56.7 

0000 2 1+ 

l.Ole-03 

3.2 

104.4 

0000 3 2 + 

1.60e-04 

10.1 

105.2 

0000 3 2- 

1.81e-04 

10.4 

115.5 

0000 3 1- 

2.88e-06 

3.9 

116.3 

0000 3 1 + 

2.61e-06 

3.9 


have assumed the rate value of the same rotation states within 
the fundamental vibration state. Collision rates values are not 
critical since, for the temperature and densities considered, T 
= 120K, 10^ cm“® < nH 2 < the transitions be¬ 

tween two different vibration states are largely dominated by 
the radiative rates. The line frequencies and the Einstein co- 
efhcier its ar e taken from the GEISA data bank (see Husson 
et al. 1992). As each of the two molecules considered, H 2 O 


and NH3, have an ortho and para configuration which are 
not coupled, we have actually run the model for four different 
molecular species. We have considered two vibration states for 
H2O: 000 (0 cm~^) and 010 (1620 cm~^), three in the case of 
NII3: 0000 (0 cm-i), 0100 (950 cm-^) and 0001 (1637 cm"^). 


4 . Results 

The population fraction of the rotation levels in the ground 
vibration state are presented in Tables a, b and ^ a, 
b for column densities of 10 ^® cm“^(A:m/s)“^ {H2O) and 
10 ^® cm“^(A:m/s)“^ {NH3), and statistical ortho to para ra¬ 
tio of 3:1 for H2O and 2:1 for NH3. Given the typical hydro¬ 
gen column density of a PDR, 2NH 2 — 10^^ cm~^(kT n/s]~ ^ 
(2 < Av < 5 over 5 km/s, see Tielens and Hollenbach 198E| ), 
such molecular column densities correspond to abundances of 
10“® and 10“® for H2O and NH3 respectively. 
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Table 5 . Calculated intensities of the brighter H 2 O lines. (T = 120K, nH2 ~ 10^ cm“®, N H2O ~ 10^^ cm~^(km/s)~^, or- 
tho/para = 3). Column 1: wavelength in micrometer; column 2: Einstein coefficient; column 3: line intensity; column 4: ratio of 
the level population fraction to the population obtained with a pure collisional model. 


A 

[urn] 

Aij 

F 

\W sr~^ (km / s)~^\ 

F 

/F(coll.) 

Upper 

V JK-K+ 

Lower 

V JK-K+ 

a) ortho — H 2 O 









538.289 

3.6e-03 

2.51e-10 

0.9 

000 

1 

1 0 

000 

1 

0 1 

273.193 

1.7e-02 

l.OOe-09 

59.2 

000 

3 

1 2 

000 

3 

0 3 

259.984 

3.9e-03 

2.63e-10 

64.7 

000 

3 

1 2 

000 

2 

2 1 

257.791 

2.4e-02 

1.30e-09 

225.2 

000 

3 

2 1 

000 

3 

1 2 

180.488 

3.2e-02 

2.21e-09 

25.9 

000 

2 

2 1 

000 

2 

1 2 

179.526 

9.7e-02 

3.37e-09 

2.3 

000 

2 

1 2 

000 

1 

0 1 

174.626 

7.4e-02 

1.89e-09 

7.3 

000 

3 

0 3 

000 

2 

1 2 

75.380 

4.8e-01 

-4.17e-09 

-119.2 

000 

3 

2 1 

000 

2 

1 2 

6.642 

7.4e+00 

4.51e-08 


010 

1 

1 0 

000 

2 

2 1 

6.185 

l.Oe+01 

-4.84e-08 


010 

1 

1 0 

000 

1 

0 1 

b) para 

- H 2 O 









398.643 

7.3e-03 

1.07e-09 

62.1 

000 

2 

1 1 

000 

2 

0 2 

303.459 

l.Oe-02 

1.87e-09 

7.4 

000 

2 

0 2 

000 

1 

1 1 

269.273 

5.7e-02 

1.50e-09 

2.1 

000 

1 

1 1 

000 

0 

0 0 

243.972 

1.9e-02 

1.68e-09 

129.5 

000 

2 

2 0 

000 

2 

1 1 

138.527 

1 .8e-01 

1.98e-10 

1.6 

000 

3 

1 3 

000 

2 

0 2 

100.983 

4.5e-01 

-2.32e-09 

-226.0 

000 

2 

2 0 

000 

1 

1 1 

67.089 

1 .8e+00 

1.81e-10 

423.5 

000 

3 

3 1 

000 

2 

2 0 

46.484 

3.0e-02 

-2.47e-10 

-2.4e+04 

000 

3 

3 1 

000 

2 

0 2 

6.672 

6.3e+00 

2.34e-08 


010 

1 

1 1 

000 

2 

2 0 

6.390 

3.0e+00 

1.16e-08 


010 

1 

1 1 

000 

2 

0 2 

6.116 

2 .1e+01 

-3.80e-08 


010 

1 

1 1 

000 

0 

0 0 
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Table 6 . Calculated intensities of the brighter NH^, lines. (T = 120K, nH 2 
ortho/para = 2). Same as Table 


10 ^ cm ^, NNH 3 


10 ^® cm ‘^{km/s) 


A 

Aij 

F 

F 

Upper 

Lower 

[urn] 


\W m~'^ sr~^ {km / s)~^\ 

/F(coll.) 

V JK± 

V JK± 

a) ortho — NH 3 





523.670 

9.5e-03 

3.15e-10 

2.2 

0000 1 0+ 

0000 0 0- 

246.773 

6 .1e-02 

6.09e-10 

4.6 

0000 2 0- 

0000 1 0+ 

169.996 

1.7e-01 

1.17e-10 

2.8 

0000 3 0+ 

0000 2 0- 

127.108 

8.4e-02 

2.92e-10 

2.6 

0000 4 3- 

0000 3 3+ 

124.648 

8 .8e-02 

2.71e-10 

3.7 

0000 4 3+ 

0000 3 3- 

11.209 

l.Oe+01 

9.09e-09 


0100 1 0+ 

0000 2 0- 

10.507 

3.1e+01 

-9.64e-09 


0100 1 0+ 

0000 0 0- 

b) para 

-NH 3 





256.575 

2 .0e-02 

1 .20e-10 

2.2 

0000 2 1- 

0000 1 1+ 

246.694 

2.3e-02 

1.25e-10 

2.2 

0000 2 1+ 

0000 1 1- 

169.990 

7.5e-02 

5.92e-ll 

3.5 

0000 3 1- 

0000 2 1+ 

169.967 

4.7e-02 

8.91e-ll 

2.7 

0000 3 2+ 

0000 2 2- 

165.728 

8 .1e-02 

5.92e-ll 

3.5 

0000 3 1+ 

0000 2 1- 

165.596 

5.0e-02 

9.08e-ll 

2.6 

0000 3 2- 

0000 2 2+ 

11.460 

4.5e+00 

6.04e-10 


0100 2 1- 

0000 3 1+ 

11.011 

5.1e+00 

6.46e-10 


0100 2 1+ 

0000 3 1- 

10.289 

8 .2e+00 

-1.03e-09 


0100 2 1- 

0000 1 1+ 

9.925 

9.2e+00 

-1.13e-09 


0100 2 1+ 

0000 1 1- 
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X (|J,m) 


Fig. 2. Synthetic dust plus H 2 O spectrum emerging from the 
PDR (dashed lines for para species). T = 120K, = 10^ cm“®, 

NH 2 O = 10^® ortho/para = 3. 


Except t he p eculiar case of the NHz metastable levels (J = 
K, see Sect. 4.2), the population fractions are far from thermal 


eqnilibrinm, because the density used does not allow collisional 
excitation rates to overwhelme spontaneous radiative emission 
(critical H 2 density ~ 10 ® cm~^ for excitation of a 100 cm~^ 
level of NII 3 or II 2 O in a lOOK gas). For this reason we com¬ 
pare the results of our model to a pure collisional model where 
the dust flux has been set to zero (Column 4 in Tables b 
and ^a, b). 

The predicted intensities of the brightest 7720 and Nils 
lines are listed in Tables ^,b and |^,b. Fig. ^ and ^ display 
the calculated synthetic infrared spectra emerging from the 
PDR. In these plots, the height of each line above or below 
the continuum is not proportional to its integrated intensity, 
because the different optical depths of the lines induce different 
line widths. The insert shows the details of the expected NIR 
vibration transitions. 

Our results show that the net effect of the NIR pumping 
is to slightly reduce the fraction of molecules in the lowest 
rotation level, increasing the population of the excited levels 
and thus increasing the intensities of the submillimeter and 
far infrared lines. However, the detailed behaviours for 7720 
and NI/s are not the same, because of the different radiative 
selection rules. 


4.L H 2 O 

H 2 O is an asymmetric top with a dipole moment of 1.94 De¬ 
bye. We refer to the rotation levels of H 2 O using the general 
numbering of the asymmetric top molecules, where J is 
the total angular momentum and K~ (respectively 77"*’) is the 
projection of the angular momentum on the symmetry axis for 
the limiting prolat e (respectively oblate) molecule (see Townes 
and Shallow 197f). The selection rules are, AJ = 0,-l-l,—1 


and A{K~ — 77+) even. The later rule forbids para/ortho ra¬ 
diative transitions since K~ — K'^ is odd for ortho and even 
for para — H 2 O. 

As can be seen from Fig. the excitation scheme for 
ortho — H 2 O is quite complex because of the combination of 
the radiative pumping by NIR photons in the vibrational tran¬ 
sitions and FIR photons in the rotational transitions. The 


cnr'l . 
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b) para H2O 


220 



vOOO 


Fig. 3. Levels and main radiative transitions for a) ortho and b) 
para H 2 O. The wavelength of each transition in micron is indicated 
in bold characters. (Thick arrows for the stronger transitions). 


molecule is mostly excited by NIR pumping in the uOlO lio 
level and subsequent deexcitation in the 221 rotation level of 
the vibration ground state (absorption at 6.18 pm and emis¬ 
sion at 6.64 pm). The population fraction of the uOOO 221 level 
is considerably increased with respect to the comparison model 
(~ 100). Two different de-excitation channels are possible from 
this level. The first one through the 108 pm( 22 i — lio) line, 
and the second one through the 180.5 pm( 22 i — 212 ) line. The 
first channel is not efficient because for the parameters we have 
selected, nH 2 ~ lO'* cm“® and ^ ~ 10 ^® cm“^(fcm/s)“^, 
emission in the 108 pm line is balanced by radiative pumping 
by FIR dust photons in the same line. Thus, the de-excitation 
occurs mainly trough the 180.5 pm line. The intensity of this 
line is increased by a factor of order 25 with respect to the 
collision model. Once in level 2 i 2 , the molecule may either be 
excited in the 821 level by absorption of a 75.4 pm dust pho- 





























ton, or de-excited in the ground level (loi) by emission at 179.5 
fj,m. The population fractions of the J=3 ladder are increased 
by one to two decades with respect to the collisional model. 
Consequently, the intensities of the FIR and submillimeter de¬ 
excitation lines are enhanced by similar factors. NIR pumping 
at 6.18 fim has thus a major effect on the (221 — 212 ) line at 
180.5 fim, whereas FIR pumping at 75.4 fim mostly enhances 
the following transitions: (821 — 812 ) at 258 /rm, (812 — 221 ) at 
260 fim, (812 — 803 ) at 278 , and (803 — 212 ) at 175 /rm. 

The intensity of the 108 j-im line is driven both by NIR and 
FIR pumping, and this line can appear in emission or absorp¬ 
tion against the dust continuum, depending on the radiation 
fielT the gas density and the molecular column density (see 
Fi^). 

The case of para — H 2 O is quite similar to ortho — H 2 O. 
The main NIR absorption line, (Oqo — lii),is located at 6.12 
am, with re-emissions predicted at 6.89 and 6.67 pm (see Fig. 
^), populating mainly the J=2 rotation ladder. Excitation 
of the J =8 ladder is much less efficient than in the case the 
ortho species because the allowed transition has a low Ein¬ 
stein coefficient (8 10“^ against 0.48 for the 75.4 pm line of 
ortho — H 2 O). NIR and FIR pumping are thus responsible for 
bright de-excitation lines in the J=2 ladder: (220 — 2ii) at 244 
pm, (2ii — 202 ) at 899 pm and ( 2 o 2 — In) at 808 pm. As for 
the ortho — H 2 O 108 pm line, the (220 — In) line at 101 pm 
will appear either in absorption or emission on the continuum, 
depending on the physical parameters of the region. 

In the regime where the NIR pumping line is not satu¬ 
rated ( Nh 20 < 10 ^® cm“^(fcm/s)“^) the absorptions at 6.18 
and 6.12 pm give a direct estimate of the ortho and para 
column densities. Furthermore, FIR lines will increase with 
the pumping rate, i.e. the radiation field. This later param¬ 
eter can be fairly well known by direct observation of the 
dust and PAH NIR spectrum. Once this is done, the inten¬ 
sities of the FIR lines can be used to trace the molecule col¬ 
umn densities. This is shown in Fig. ^ where the calculated 
iso-intensities of the 180.5 pm (221 — 212 ) line are drawn ver¬ 
sus the H 2 O column density and the H 2 density. The dot¬ 
ted contours in the same plot correspond to the comparison 
model (collisions only). In the comparison model the line in¬ 
tensity increases fairly regularly with the gas density and the 
molecule column density. If radiative pumping is introduced as 
in our model, it dominates the excitation at moderate densi¬ 
ties, n_tf 2 < 10 ® cm“®, so long as the 6.18 pm pumping line is 
optically thin: N h^o < 10^® cm~^{km/ s)~^ . In this domain of 
the parameter space, which is appropriate for the conditions 
that prevail within PDRs, the line intensity is almost indepen¬ 
dant of the gas density and can be used to trace the column 
density. However, for the 180.5 pm line, the intensity is not a 
linear function of the column density. This is because the upper 
and lower levels of the transition have several other radiative 
exciting or deexciting channels which become optically thick 
in different column density ranges. 

Fig. ^ shows the intensity contours of the (221 — lio) 108 
pm ortho — H 2 O line. For the sake of clarity we didn’t over¬ 
plot the contours obtained with the comparison model. They 
are quite similar to those of the 180.5 pm line in Fig. 

As mentionned above, this line appears either in absorption 
or emission, depending on the gas density, the molecule col¬ 
umn density and the radiation field. Including NIR pumping 
favors emission in this line since the upper level is directly 
fed by the de-excitation of the excited vibration state. This 



X (|J,m) 

Fig. 4. Synthetic dust plus NH 3 spectrum emerging from 
the PDR. Same as Fig. 0 . T = 120K, = 10^ cm“®, 

htNH 3 = 10^® cm“^(A;m/s)“^, ortho/para = 2. 


is why the domain where the line will appear in absorption 
against the dust continuum is restricted to a small region: 
10^“^ < NH 2 O < 10^® cm“^(A:m/s)“^, lO’^ < Nh 2 < 10® cm“®. 

It is worth pointing out that, in addition to increasing con¬ 
siderably the H 2 O submillimeter and FIR line intensities, ra¬ 
diative pumping will strongly modify the lines intensity ratios 
with respect to pure collisional models. For instance, the 180.5 
pm to 179.5 pm ortho — H 2 O line intensity ratio is increased 
by a factor of order 10, and the 898 pm para — H 2 O to 588 
pm ortho — H 2 O submillimeter lines ratio is enhanced by a 
factor of 60. Such effects have very important observational 
consequences, particularly when the lines are used to derive 
column densities and ortho to para abundance ratios. 

Finally, the NIR lines predicted to appear in absorption and 
emission against the PAH continuum can directly be used to 
probe the H 2 O column density. The line to continuum contrast 
of the emission transitions is high (see insert in Fig. because 
the intensity of these lines comes from energy absorbed at the 
maximum of the 6.2 pm PAH band. With ortho — H 2 O, we 
actually have a case of quasi-resonant excitation. The primary 
emission transition at 6.64 pm is thus a good target for ob¬ 
servations. For instance, if we assume a line width of 10 km/s, 
the measured line to continuum contrast is still of order 50% in 
the 6.64 pm transition, after dilution by the ISO-SWS grating 
spectral resolution {R ~ 1500). This is to be compared to the 
contrast in the absorption line at 6.18 pm, ~ 5% for the same 
resolution. 


4.2. NH3 


NH 3 is a symmetric top prolate rotator with a dipole of 1.47 
Debye. The possible tunelling of the nitrogen nucleus through 
the hydrogen plan is responsible for the inversion splitting of 
the rotation levels with K ^ 0 (levels + and -). We refer to the 
leve ls usi ng the usual J,K± numbering (see Townes and Shal¬ 
low 1975). An additional quantum number due to a rotation- 


vibration interaction, I = ±1, appears in the case of the degen¬ 
erate vibration states (e.g. 0001). The selection rules are those 
of a symmetric top: AJ = 0, -|-1, —1; -f ; and AK = 0 for 
transitions which do not involve a degenerate vibration state. 
The later rule implies that the different K ladders are not ra- 
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cm'l 

1050 


950 


100 


0 



by a factor ranging from 2 to 5: (4,3 —3,3) at 125 and 127 fim, 
(3, 0 — 2, 0) at 170 /rm, (2, 0 — 1,0) at 247 nm and (1,0 — 0, 0) 
at 524 pm for ortho — NH 3 ; (3, 2 — 2, 2), (3, 1 — 2,1) at 166 
and 170 pm and (2,1 — 1,1) at 247 and 257 pm concerning 
para — NH 3 . 

As for H 2 O, the effect of NIR pumping is relatively more 
efficient for the transitions from the higher rotation levels of 
the downward cascade, because it is more difficult to populate 
these levels by collisions with the gas. However, in NH 3 the 
metastable levels are thermally populated whatever is their 
energy. This facilitates the collisional excitation of high energy 
levels, if they are not far from the bottom of the K ladder . For 
instance, levels (4,3)± at 165 cm~^ are only 79 cm~^ above 
the (3,3)± levels from which they can be efficiently populated 
by collisions. This explains why NIR pumping is relatively less 
efficient for NH 3 than H 2 O. 

Finally, there is no complex behaviour in the NH 3 lines as 
there is in H 2 O (e.g. the 108 pm ortho — H 2 O line). This is 
because each K ladder is the location of a unique de-excitation 
cascade, forbidding competition between different channels as 
in the H 2 O asymetric rotator. 
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vOOOO 
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Fig. 5. Levels and main radiative transitions for ortho and para 
NH 3 . Same as Fig. 


diatively coupled. Consequently, the lower -|- and - rotation 
levels of each K ladder (J = K) are the bottom of a radiative 
cascade, and are almost in thermal equilibrium with the bot¬ 
tom levels of the other ladders. They are usually refered to as 
metastable levels. 

Three vibration states have been included in the excitation 
calculation of NH 3 : 0000, 0100, and 0001. However, despite 
comparable Einstein coefficients, the frequency dependency of 
the radiative pumping (see formula ^ favors radiative pumping 
in the 0100 state (950 cm“^) with respect to the 0001 state 
(1600cm“^). The 0100 vibration at 950cm”*^ correponds to 
the oscillation of the nitrogen nucleus above (or bellow) the 
hydrogens plan. 

The NIR pumping induces a pump cycle within each K 
ladder, with several NIR absorption and emission transitions 
predicted between 9.5 and 12 pm (see Table ^). The intensities 
of the FIR and submillimeter deexcitation fines are increased 


5 . Conclusion 

The main conclusions that can be drawn from our study are 
the following: 

1) NIR photons are able to contribute very significantly to 
the rotational excitation of polar molecules within PDRs in 
galactic star-forming regions and starburst galaxies. 

2) As this mechanism is independent of the gas density, 
we predict that bright FIR and submillimeter rotation fines 
should be observed from the relatively moderate density and 
non shocked gas in PDRs, nH 2 < 10®cm“®, where molecules 
are expected to be abundant because of grain mantles evapo¬ 
ration. 

3) For H 2 O and NH 3 the intensity of the FIR and sub¬ 
millimeter rotation fines can be increased by a factor ranging 
from 2 to several 100 with respect to models which do not take 
into account NIR pumping. This means that the molecular 
column densities that can be infered from observations toward 
such regions, using standard models, can be wrong by similar 
factors. 

4) For gas densities below 10® cm“®, where NIR pumping 
dominates over collisional excitation, the intensities of the FIR 
lines becomes mostly independant of the gas density. This im¬ 
plies that the FIR and submillimeter lines can be used to probe 
the molecule column density if the radiation field has been cor¬ 
rectly estimated. 

5) For ortho — H 2 O we have a quasi-resonant excitation 
of the molecule, by absorption of 6.2 pm PAH photons in the 
6.18 pm fundamental vibration transition. We thus predict a 
high fine to continuum contrast for the emission de-excitation 
line at 6.64 pm. 


NIR pumping is likely to affect not only the PDRs, as mod¬ 
elled in this paper, but also the molecular clouds and cores in 
their vicinity. The extinction in the 6 to 8 pm range is mini¬ 
mum (typically 0.02 of visible extinction, Rieke and Lebofsky 
and thus, the radiation produced in the PDR penetrates 


198! 


well in the molecular clouds. 

Most molecules present in PDRs and nearby molecular 
clouds will actually be affected by radiative NIR pumping. 
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a) (221-212) 180.5 |im, 0-H2O 



b) (221-1 10 ) 108 |im, 0 -H 2 O 



Fig. 6. iso-intensities of a) ortho — H 2 O {221 — 212 ) 180.5 iim and b) ortho — H 2 O (221 — lio) 108 fim lines as a function of the molecule 
column density and the hydrogen density. Dotted lines in a) are the contours obtained with a pure collisional model. Dotted lines in b) for 
an absorption line. Units are Wm~‘^sr~^(km/s)~^. 


In addition to H 2 O and NH 3 one can list CS, HCN, CO 2 , 
CHi, H 3 O+, HCsN, and CH 3 OH. The combination of the ab¬ 
sorption and emission lines from all species existing in the gas 
phase is likely to be responsible for the very complex detailed 
structure of the 5 to 20 /rm spectra measured at moderate 
spectral resolution in the direction of the M17-S W PDR with 
the ISO-SWS intrument (Verstraete et al. 1996). Decrypting 
such spectra in terms of molecular abundances will require the 
production of synthetic model spectra including a significant 
number of molecules of astrophysical interest. This will be done 
in a forthcoming paper. 
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